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PREFACE 


This paper vresents the results of a study of imposing a three- 
dimensional perturbation on a barotropic westerly current. The 
purpose of the present investigation was to present theoretical 
solutions to the following basic problems: (1) the determination 
of the three dimensional structure of the — (2) the determination 
of the speed of propagation of the wave; (3) the value of considering 
vertical motion in the wave motion, 

Underteken 2s the thesis reguirement for the degree of Master of 
Science in Aerology, this »aper wes prepared at the United States Naval 
Postgraduate School, Monterey, California during academic year 1951-1952. 

The author is particularly indebted to Associate Professor F, L, 
Martin of the Department of Aerology, for his advice and guidance during 


the entire preparation of this study. 
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I. INTRODUCTION 


E 


The theoretical study of the behavior of large-scale flow patterns 
in the extratropical zones of the earth, associated with migratory long 
waves in the belt of orevailing westerlies, occupies an important place 
in meteorological literature. 

Many informative theoretical and synoptic studies of the behavior 
of these large-scale flow patterns have been made since V. Bjerknes [5]. 
in 1916, advanced a theory, based on general hydrodynamic considerations, 
that cyclones originate as dynamically unstable wavelike disturbances in 
the westerly current. In 1923, H. H. Clayton [27] made the first 
systematic use of zonal index by means of the trend of mean pressure 
versus latitude to determine increasing and decreasing zonal circulation. 
By 1933, V. Bjerknes and colleborators L27) had computed the zonal 
velocity distribution from zonal temoerature and pressure distribution 
at the ground, showing the maximum horizontal velocity normally occurs 
just below the tropovause and concluded that there were 2-3 waves around 
the 60° latitude circle and slightly above 4 around the 30° circle. In 
.ل و1937‎ 5 [5,28] offered a simple explanation for the displace- 
ment of perturbations superimposed on zonal distribution. Variations 
in latitude and pressure perturbations, with simultaneous transport of 
air, result in changes of amplitude or longitudinal displacement of 
long waves. This result was obtained by combining the equations of 


motion with the equation of continuity. Rossby Gs] pointed out two 


(1) 








basic objections, first, the omission of acceleration terms which 

may be of the same order of magnitude as the centrifugal and de- 
flecting force, and second, the displacement of observed flow patterns 
caused by isallobaric systems resulting from divergence associated with 
initial wind distribution. Rossby [27] extended the J. Bjerknes paper 
reasoning that since the depth of the atmosphere was much smaller than 
the horizontal scale of large disturbances, it was permissible to treat 
the atmospheric motion as two-dimensional, non-divergent flow. Further, 
he deduced that between the longer westerly moving waves, there existed 
an intermediate wave length, such that its perturbation remained 
stationary. The speed of propagation of moving weaves is dependent on 

the strength of the westerlies and the wave length, and was based on 

the consideration of a change in vorticity in a verticel air column 
which has been displaced from one latitude to another. Haurwitz (16, 17] 
extended Rossby's 1939 solution by studying the effects of lateral limits 
of the wave and a variable Coriolis parameter on a curved earth. 

Wamias and Clapp [22] investigated Rossby's wave speed formula 
statistically. Their verification was generally successful and they 
concluded that amplitude and zonal velocity were both important in 
determining the stationary wave length. 

V. Bjerknes had, by this tine, found the number of waves around the 
145% latitude circle to be between four and six. In 1944, J. Bjerknes 
and Holmboe (19) derived a formula enalegous to that of Rossby for the 
more general barotropic atmosphere which agreed qualitatively with that 


of J. Bjerknes! of 1937. 


(2) 
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By 1947, Gherney (5) in his baroclinic atmosphere theory, based 
on J. Bjerknes and Holmboe's theory of weve motion in a baroclinic 
atmosphere, introduced a solution to the barotropic wave. ith this 
theory, he obtained a solution, based on the principle that a wave 
will travel with such a speed that the pressure tendencies arise from 
the displacement of the pressure pattern in accordance with the field 
of horizontal divergence. Charney found this physical explanation 
acceptable, but the dynamics of the wave could not be solved by 
analytic methods based on semi-empirical considerations of gradient 
wind. Charney (1947) integrated the eouations of motion, assuming 
motion to be adiabatic, and thus determined the speed of propagation, 
stability criteria end the three-dimensional wave structure. In 1948, 
Eliassen 27] generalized the Rossby~Holmboe wave formula, but he 
neglected vertical acceleration and the Coriolis parameter when it 
appeared in multiplicative combination with the vertical velocity. 

In 1949, Cherney and Eliassen (71 extended numerical analysis by use 
of the quasi-geostrophic approximation. 

In a field where horizontal divergence was so small that the 
absolute vorticity of air parcels was conserved and the scale of 
motions was large enough that the Coriolis parameter variation with 
latitude was important in determining wave dimension, Cressnan L9J 
scrutinized the Rossby wave velocity formula and found retrogression 
in a series of quasi-stationary long waves began when the basic zonal 
velocity decreased in speed or shifted south. Such retrogression in- 


creased the long wave number (i.e., formed a new cold troush). 


(3) 










jesi md vmm m "m 
— دوم‎ ve 





E 
E a 
٧ 





Continuing his previous investigation, “hzrney iZ; 8] found that 
the mean flow with respect to pressure could be described as two- 
dimensional non-divergent in an eouivalent barotropic atmosphere, 

He used this method as a basis for partially successful twenty-four 
hour forecasts of the 500 mb flow pattern by numerical integration 
of the non-linear vorticity equation. The simplified equations 
derived expressed conservation of the vertical component of absolute 
potential vorticity. The hydrostatic equation together with the 
individual rate of change of vorticity evaluated geostrophically, 
reduced to a single partial differential eouation for pressure 
tendency. The basis of correspondence between the barotropic and 
baroclinic atmosphere lay in the notion of an equivalent barotropic 
atmosphere, in which horizontal motion at a particular level approx- 
imated the actual motion. 

It cen readily be seen that a great amount of effort kas been 
placed on the long wave theory in an effort to gain a forecasting 
tool. The use of vertical motions has suffered in comparison, but 
primarily because of the sparsity of observations and the necessity 
to rely upon indirect information. Investigations of the vertical 
component of motion were undertaken as far back as 1911, when V. 
Bjerknes [13] described a method for such computations. In 1913, 
Sxner (13) suzgested the high correlation between 9 km pressure and 


temperature could be accounted for by vertical motions. Hesselberg, 


(4) 








in 1915, extended Zxner!s argument by showing thet the oresence of 
persistent vertical motions of 1-10 centimeters „er second would 
account for the observed magnitude of pressure-temperature changes, 
Jeffreys [24] in his 1922 stuày of vertical velocities necessary to 
maintain hailstones, found vertical velocities of any magnitude at 

all occurred only in eddies and could be ignored in treating mean 
motion. In 1927, Haurwitz 13] showed that in an atmosphere of 
constent lepse rete, that the maximum vertical component of velocity 
wes in the region of méeximim divergence and was of the order of eight 
centimeters per second. Brunt, [4] in his Memoirs (1934), found vertical 
velocities involving divergence were of the oräer of a few centimeters 
per second, 

As observational information increased, quantitative analysis of 
divergence and vertical currents became more practicable. wurst (1) ٩ 
in 1940, found vertical velocities of a few centimeters per second; 
Fleagle [2], in 1945, Gaviola and Fuertes [i5], in 1947, assumed 
adiabatic processes and computed the vertical component of velocity. 
Panofsky (231, in 1946, from third degree volynomials, found values 
of the same magnitude (i.e., a few centimeters per second). 

4hile vertical motions of several meters ver second have been 
observed, their duration was short and smoothed out over a twelve- 
hour period. Thus, the high valued measurements by theodolites as 
observed by Suring [15] have not altered the basic contention that 


persistent vertical motions ere of e much smaller order, 
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Sharney 16), in assigning magnitudes to vertical motions in 1948, 
explained that the failure of the tendency equation in pressure prog- 
nostication was due to the discounting of vertical notion. Charney ts), 
further determined that forecasts with the quasi-geostrophic assumption 
could be achieved only when the distribution of vertical motion and sea~ 
level pressure distribution were know. Thus, it is necessary to extend 
numerical forecasts to permit forecasting vertical motion. In 1949, 
Bellamy [2) described a method of objective calculation of vertical 
velocities. In 1951, Phillips [25] indicated that the equivalent 
barotropic atmosphere was incapable of oroviding distribution of vertical 
motion and sea-level pressure adequately for forecasting the "weather", 

In view of the interest in vertical velocity distribution, it is the 
purpəse of the present investigation to »resent a theoretical solution 
of a two-lsyer barotropic, incompressible atmosphere with three-dimensional 


velocity divergence zero, introducing vertical motion as 2 perturbations 


(6) 
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II. FORMULATION ASD RESULTS GF YAZ PROBLEM 


å brief summery of the principal subject matter of the investigation 
unhempered by mathematical detail is now given to set forth the procedure 
followed anâ results obtained. 

Chapter III concerns the construction of a barotropic model. The 
troposphere is characterized by a constant vertical lapse rate; the 
stratosphere is assumed to be isothermal, with constant zonal wind 
throughout the atmosphere. The equations of motion are given in ordinery 
cartesian coordinates for an incompressible atmosphere. Horizontel 
velocity divergence is assumed to be compensated by vertical motion, 
Boundary conditions are formulated. It is shown that the model is 
consistent with the boundary conditions and the steady state notions 
prescribed, 

The last section of Chapter III, treating the actual motion asa 
small perturbation superimposed on an undisturbed constant zonal current, 
is now discussed. The velocity components, as well as the pressure 
component, are expressed as — harmonic perturbations of infinite 
lateral extent traveling in a west-east direction at constant velocity. 
These perturbation components are functions of height as well as of 
horizontal coordinates and time. The vertical velocity, involving a 
general expression of the z-cooräinate, is used as a basis for the 
determination of all other velocity components. To simplify the 


mathematical difficulties, the meridional perturbation amplitude is 


(7) 








made proportional to that of the zonal perturbation. The zonal 
amplitude, in turn, by the prerequisite of zero of velocity divergence, 
is simply related to the vertical motion amplitude function, H(z). It 
is shown that pressure is also a function of the vertical motion ampli- 
tude function. 

Since the vertical amplitude is in terms of H(z), the amplitudes of 
the other components are determinable from H(z). The determination of 
the proportionality factor between the meridional and latitudinal per- 
turbation velocities assumes an important role in the solution of the 
problem. Upon its determination rests the dependency of the wave 
velocity and the maximum allowable meridional perturbation velocity 
for a particular maximum zonal perturbation velocity. By analogy to 
Rossby's wave velocity, a restriction thet the proportionality factor 
be less than unity is assumed. It is shown thet the proportionality 
factor is dependent upon the discriminant of an auxiliary quadratic 
equation. Following the boundary conditions this discriminant must 
have a value slightly in excess of unity. The discriminant value 
must be determined from an inequality. The choice of the value of 
this inequality plays a dominant role in any further solution of 
this problem. 

The problem then becomes one of determining the speed of propa- 
gation of the wave and the character of the vertical perturbation 
which introduces the amplitude as a function of height. The deter- 


mination of these are dependent upon wave length and the parameters 


(9) 





characterizing the mean state of the atmosphere, namely, the height 
of the tropopause, the vertical velocity at tne tropopause, the 
zonal wind, the vertical la»se rate, the mean latitude and surface 
density, pressure, and temperature, 

In Chapter IV, the form of the perturbation is set forth, The 
perturbation functions ere solved in the stratosphere, ano. the tropo- 
sphere (by approximation). "ne proportionality constant between 
meridional and zonal perturbation velocities is determined., This 
constant is of zreat importance in the determination of the wave 
speed and stationary wave length. Amplitude evaluation is determined. 
The complete solution of the approximate case is shown diagramatically. 

In Chapter Y, the tropospheric perturbation functions are solved 
without simplifying the coefficients of the differential eouation 
satisfied by the amplitude function, H(z). 

Due to barotropy, no phase change is introduced into the structure 
of the wave. Vertical cells do &poear in the perturbation velocity 
structure moving with the basic wave velocity. The perturbation 
fields of velocity diminish with increasing height and eventually 
approach zero. The wave in the meridional velocity field legs 90° 
behind the weve in the oressure field and the wave in the density field 
is in phase with the pressure wave. The zonal perturbation is in phase 
with the oressure wave, whereas the vertical velocity perturbation leads 


the pressure wave by 90%, These phase relationships are shown in 


Figure 3. 


(9) 
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III. Toe MATCEESMATIOxL PORMOLETIOX 


l. The Atmospheric Model. 

AS an approximation to the atmosphere in middle latitudes, e 
model of the undisturbed state is characterized as follows; (a) the 
undisturbed wind is zonal; the speed is constant and independent of 
the coordinates throughout the troposphere ənd stratosphere; (b) the 
lapse rate of temperature is constant in the troposphere and zero in 
she stratosphere; (c) both troposphere and stratosphere are barotropic, 
incompressible and frictionless layers. 

2. Fundamental Equations. 

It is assumed for purposes of mathematicel simplicity that the 
curvature of the earth cen be neglected. This assumption is valid when 
length of wave is small compared with circumference of the zonal circle 
along which the wave moves. ‘This planar motion is expressed ina 
rectangular system of coordinates x, y and z with x increasing to the 
east, y northward, and z vertically upward. Then u, V, W denote 
corresponding total velocity components (consisting of undisturbed 
velocity and perturbation velocity). Further, denote totel pressure 
by D, anguler velocity of the earth byl, acceleration of gravity 
(assumed constant) as g, geographical latitude by 9 , the z and y 
components of Coriolis acceleration, f and 73 respectively. The 


aensity JO , having been assumed incompressible, is constant and 


(10) 





equal to the Unperturbed value /4 .he sulerian equations of motion 


in ordinary cartesian coordinetes, fixed in the earth, become 


a. 85 ፦ ፖ - “3.2 (1) 
رال ووي‎ B c 

dž — Lab - fu 

dt A dY ) (2) 
dur - MRE gü — 

ur “Pz Z 3 * (3) 


A fourth equation is obtained from the assumption that the three- 


dimensional velocity divergence is zero. 


Y, Y 20 , (^) 


A relationship between density and height in the troposphere (with 
the assumption that g is constant with height) provides the fifth 
eguation 

I 


yo ور‎ AE (5) 


But the density-height relationship in the stratosphere is 


E (2-hr) 
2 ዳጊ (6) 


A 


ሥሥ = #5 aT 


where An is height of tropopause. 


(11) 
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3. The 2oundary Conditions. 
The bou- dary conditions express the following ohysicel oroperties 
of the motion: (a) the perturbation components are continuous ecross 


the tropopause, that is, 


AB ME AV = Al - Ade AT 6۵ (7) 


(b) the momentum vanishes at the limit of the atmosphere. 


lim fe = him. መሠ = lim pus z NOME (8) 


230 73 مه‎ 
4, Derivation of the Perturbation Zouation. 

The fundamental eouations (1-6) together with the boundary con- 
ditions (7), (8) impose the necessary restrictions on the theoretical 
model.  ;ith the convention of a capital letter for steady state value 
and a small letter for the same quantity in the perturbed state, the 
characterization of the steady state is given by the conditions 


V = qu = OP. - © , U = eomsbant. (9) 
WW: 2r OX 


Equation (2) implies thet the undisturbed flow must satisfy the 


condition of geostrophic eouilibriun 


- AME المت‎ (10) 
24 


Souation (10) implies that jo and P sre functions of y and z only. 
Xouetion (3) exoresses the condition for hydrostatic equilibrium in 


tne mean state 
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with the omission of BU, which can be shown both empirically end 
theoretically to be negligible in comparison to the acceleration of 
grevity. 

By means of the equation of state, end equations (4) and (5) with 
prescribed values of U, Py Ty, FT, and hp, the theoretical steady state 
model is completely described. "he boundary conditions(7) end (8) are 
also satisfied by our steady state model. Hence the model prescribed 
is consistent with the fundamental eovations (1-6) and the boundary 
coadi tions (7), (8). 

5. The Linearized Houations. 

The motion to be investigated can be treated as a small perturbation 
with velocity components u, v, and w superimposed on an undisturbed 
constant zonsl current U, with the assumption of infinite leteral 


extent. ‘Thus the motion is described by 
U = u(x,z,t)+ U, V = v(x,2,t), 9 = w(z,z,t) 5 )12( 


Similar expressions for the density end pressure in the disturbed state 


are 
"mz بيهو ) م5‎ 6 ( + END, P= ሥ(።›። > (13) 


In deriving the perturbation equetions the usuel assumotions are 
made. (1) The total disturbed plus undisturbed motion satisfies the 
hydrodynamic equations, 2s well es the undisturbed motion alone. 


(2) The verturbations are so small that terms of second order in the 


(13) 
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perturbation quantities can be nezlected with resvect to terms of 
the first order in the >erturbstion quantities.  Lineerizing the 


Sulerian equations of motion (1-3) and simplifying by 


ره 


meens of the 


steady state equations (9-11), the system (1-3) becomes: 


Ju + او لا‎ + Bw - fvz-! Op (14) 
ot 2x هر‎ 2 


av + Vou +fu = O (15) 


DE Ox 
0 عن‎ Udur — wù :-L2 (16) 
ጋቲ M Ox ? 2 $F 


du TE = ©, )17( 


(14) 
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IV. THE BAROTROPIC PERTURBATION 


l. The Form of the Perturbation and Resulting Conditions. 

The verturbation components, traveling in the x- direction at a 
constant velocity, are assumed to be simple harmonic motions of infinite 
lateral extent. The amplitudes of the perturbationare included in a 
general expression of the z-coordinate, as follows 

u = sin m(x-ct) F(z); vT = cos m(x-ct) Q(z); 


(18) 
= -m cos m(x-ct) H(z); p 7 sin m(x-ct) G(z). 


= 
i 


To simplify the mathematical difficulties, the assumption is 
made that F(z) is proportional to Q(z), D being the constant of 
proportionality. Since u and w of eouations (18) are related by the 


condition of zero velocity divergence, it can be shown that 
F(z) = ፪! (2) = Q(z) D. (19) 


The equation for the velocity of the wave is derived by substituting 
the perturbation velocities (18) into the equation (15). #ith the 


equalities of (19) known, the wave velocity equation is 
c = U - fD (20) 
m 


Since U and m are real, c will be either real or complex, depending 


on D. The determination of 2 will be discussed fully in Chapter IV-2. 


(15) 
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Substituting the perturbed velocity, pressure and mess fields into 


the perturbation eguations (16) end (14) leads to 


bai ef Dino Zä Hai. (21) 


GE) مر‎ SH H(z) t23 H'tz), (22) 
Dm 
respectively. 


The differentiation of (22) with respect to z, leads to 


/ 
zl- DŁ 7 ተ(ጅ፡- / ፉን z). (23) 
Gz) [| H(z) cl e Hat pBH tte 4 Ha) 


By differentiation of density from equation (5), we have 


“/ 2 
ge = TG a) : (24) 


Equating (21) and (23) and replacing /ወ end g2. by their 
2 
expressions in (5) and (24) leads to the following ordinary 


differential equation for the amplitude function H(z): 


- «z) Ht -KHO + - p*ri(i-«2) H@)=9 (25) 
1 )He) -K ën Cpt) ODED : 


where K is an abridged notetion for 


ff=< (3ARr-)， (26) 


Using the further abridged notations, 


Az Akpm B= — Dim” (27) 





(16) 
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equation (25), the troposoheric differential equation becomes 
ሪ ( gá-«2] He 7o - 
፦-=2) ዘሂጋ-- ደ ዘር3 ተኒሳተ (28) 


Ze Perturbation Functions in the Stratosphere, 

In exactly the same manner as described in the —À section, 
the stratospheric differential equation follows by using (6) instead 
of (5) for the density expression. Thus, analogously to equation (24), 


the stratospheric density derivative is ١ 


== ES - p/ h, ) (29) 


where hg = RT,/g is the height of a homogeneous layer at stratospheric 
temoerature 22 4e, therefore, obtain the second order linear differ- 


ential equation with constant coefficients 


2) — H! o. 
HG io ar —5 Hé = (30) 


The auxiliary eovation of (30) is 


(31) 


— 了 i 
whence 


I-42Dmh. ز20‎ 
aep =) 


Mhe determination of the value of D is dependent upon the value 


(22) 





essigned to the discriminant of the auxiliary equation (32). By 


analogy with the Rossby wave velocity 
= p 2 
U « «A/Em 


(17) 





"ei 


it would follow that, for Rossby waves, 


Dp = /3/Efm Z o. 


The only restriction placed on D at this stage of the investigation 5 


is that DS 1, from which it follows at once, that the discriminant of 
(28) has a velue n at latitude 459, slightly in excess of unity, 


n,» ገቺ 


Solving the discriminent for D: 


DL (emh Y- 1) |- 0 +(n2-1) =0, (33) 


let 
2 = (27 -1) Y E (amh,) | (25) 


Solving equation (33) by the quadratic formula gives 


& መሼ —— SEA, P 


D 7 (n?-1) 34)ኙ Ta ay 


In order to sxpand the radical in (25) by means of the binomial 


expansion, the condition that 


2 
O8‏ دزم مھ د( ځ)/ للا » | 1- Re‏ 


must hold. 


* D Z1 with a moderately long wave of 6000 km, would give & value 
of U-c >» 100 mos, so that obviously D 21 does not apoly. 


(18) 








INEO EE b A 


we OSS ፥ 
emmi <” NN came lade '* 4 p a 
A ጮጮ Wm Om Ap ir, $ 


به 


ue tmt اد گرم ټس‎ ይ 





ee ee 








غلا 
du Tum Ya a "5‏ 
ዎነ አህ”. .፣= * "ው Art‏ 06 
1 
Tl IET. o, iğ‏ - 
m Cada ١ Y | di‏ 


(361) 


(37) 


(28) 


Zxamining the two possible cases: (a) GEE ODO ETT 
Case (a): Y, <1. 


mquation (36) becomes 


(2 mn >(ክ:-13./' , 


so that 


2. 
For L = 6000 km, n4 < 1.012. Then from the inequality of 
case (a), the expansion of the radical of (35)with the choice of 


the negative sign, yields 


D = (n-I = 十 十 DY. 
AX š E" 


Oase (b): Y,» 1 


Equation (34) becomes 


(2m h) /(r£-1) 21. 


It follows thet 
> ل‎ 
| < ረ Cis nay 2 


where n „< 1.0004 for a wave length of 6000 km. The value of mec 


determined by case (b) is included in that of case (a), so that >يلا‎ 1 


will be used in the evaluation of D, Moreover, all terms of (37) beyond 


the second will be considered negligible. Thus at latitude 45° where 


رک بت کر 


(12) 





$ 


(39) 


„Ch mh | 


qid 


D=4(m-ı E 
2m hz 





The evaluation of D, which for a particular mean latitude and wave 


number remains constant, is dependent on the value of n „selected 


In this investigation, it was assumed 


from the inequality (36!). 


desirable to have the ratio of the absolute maximum velocity of the 


horizontal perturbation to that of the vertical perturbation as large 


(29!) 


The value of D is com- 


This ratio, hereafter called 


as consistently possible. 


M = Home| = $, (I~ nà 
yn 


UJ max 


can be shown to be à function of n , and L, 


pared in Table 1. for various values of n, and wave number i. 


376 
238 
wre 
.377 
.377 
. 377 


0.5 
0.5 
0.5 
0.5 
0.5 
0.5 


= 0 
y + 218% A 


m in 
meters ~ 
N Ne x 10 2m ha — 
4 1.0056 0.888 0.0113 05107 
5 1.0070 15210 0.0141 0.119 
6 1.0084 127392 0.0170 0.132 
7 1.0098 1.554 0,0198 0,141 
8 120112 1.776 0.026 0 
9 140125 1.998 0.0254 0.161 


Table 1. 


(20) 


Constant parameters: Ø = 450, T 





Having determined the value of D qualitatively in terms of m, 


h, and n the general solution of the stratospheric equation (30) is 


22 


He) > e መጫሬ — $22 E Cle + Sz nZ : (አዕ) 


with the abridged notation 
S, 一 V2 ho ። 
By application of the boundary condition (5), our -?0, it follows 
27% 
that Haro if and only if the coefficient C, is made zero. 
Z-»*o 2 
The particular solution satisfyinz the boundary conditions is 


Hz) > E ۱ 


yy 


In the determination of the pressure perturbation in the stratosohere, 
the amplitude function G(z) is rewritten from eguations(22) and (6), in 


the abridged notation, as 


zl 
ሀሪ) = G 2 (SFU-DIU-n] +3) (1- «d, KT CR, 
fo G — ال‎ Yi ሣን Lë 


From the solution of (41), the perturbation components of the 


stratosphere are defined. In abridged notation, they are 


um ር (,م-/) ري‎ e lte SIA m c-cz), 


$, 0-4 


VI CS OE Cos m(x-c&), (5) 
D 


(ተክ) 
Ww, = - 0 C, e* 2) Gas MCX-cH, 


p, = Ç G) sin my-cb. 


(21) 
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In order to solve (43) for Cj, suitable numerical values must be 
assigned to the constant parameters €, hg, T, + We, rand n; AE. 
The following values are selected: Y = 499" مط‎ = 10 km, T, “ 219656) 
Wmax = 5.5 om sec^l, 8 = 5, 83 = 1,007, کش‎ = 0.001121 gm em. 
With assignment of these quantities, numerical values may be derived 
for the parameters m, S5, M, where m = 1.11 x 1076 meters^l, 
53. = ,0785 kml, M = à. Zquation (43) yields Cy 7 -5 x 10“ m“ sec”2, 
3, An Approximate Solution in the Troposphere. 

To obtain the solution of the barotropic wave by approximating 


the value of (1 -«z) in the tropospheric equation, (28) let 
Caz) رر‎ (Bh) 


be a satisfactory approximation. Substituting this approximation 
into equation (25), the equation becomes a second order differential 


equation with constant coefficients, namely, 


Spe A tt Cz) ተ (Are) Hcz) = 0. (45) 


o 


— 


ር) 9/0 | 06 . | (46) 


Then, if h, = RT; lo ‚ equation (45) can then be shown to be 





Hay) -&_ Hm +(BEDm_ _ SS: ga) «o. Un 
Hs e h..{p*- |] —— 


(22) 
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The general solution of this equation is dependent upon the discriminant 
obtained from the auxiliery equation of equation (47). By letting the 


discriminant equal h, » where 


2 
= | — 4/3 h; F 4 (0 m h,) ) (43) 
0 FED) ሮሮ l) 


and introducing the abridged notation 53 = € /2. h; ) 
with D known from the stratospheric solution, 3 is determined. 


Then the general solution in the approximate case is 


H(z) = > [C, e57 + ee | | (49) 


h, Application of the Tropovause Boundary Condition. 

The method of evaluating the coefficients C, and Cg is the same 
for both the approximate and actual case. The approximate case is now 
considered. 

with the boundary condition of a discontinuity of the first order 
at the tropopause the perturbation components must be continuous across 


the tropopause, 


UL, 2 03 P WU, 0X ) — بت تت‎ )50( 


Then from (49) the first derivative with respect to z is 


— SE 
NO - Ser [Gs n9 e^ "+ (‹,( “ክነ ያሚ s] (51) 


(22) 








so that 


-n,9 
ን ከርሮ h; >) Coş m-ct), 


uz” S,e5s*[ C, Gras) e "2%* 4 ‹፡(/-ዛ3)67 55] 547 m (-c£), 
n 2 -nz S. 
ma)” T me IC Ce 7 +C, e : e 


5,2 
U 3C max) z 93 € ° [ein E Ca Ci- naye Be | - 


me (Ge E‏ ~= وس 


(52) 


From the boundary condition (50), equate vertical velocities and also 


horizontal velocities at the tropopause, thus 


C e S. Ci- nah > Sah? (( eS, hr #C,e-%55s hr) 
> 


(53) 
C, 9, (“ነ esk na) hr. Ses sale. nee, Ca -n;) etie] 


By solving these two equations simultaneously Ca and C, are determined. 


3 
5, The Complete Solution in the Anproximate Case. 


With C, and Cz determined, the »erturbations of the troposphere 


ere by approximation defined as 


43 =  Ugimax) Sin m(x-<t), 


` 


D “3(maz) “es m-et), 
(54) 
uy =< “3 Cmax) Cos mCx-ct), 


e = ዕ- «27 ' L Be U3 (max) - ممح وس‎ | Sın mc L), 


ESCH 





In order to solve (55), suitable numerical values must be 
assigned to the parameters of the troposohere, in addition to those 
assigned in the stratosphere (e.g., Section 2 of this Chapter). The 
ađditional parameters, having selected values, are 
Tg = 256,59A,[* 2 5,180 kml, with assignment of these quentities it 
can be shown thet the constants C4 and C, have the values 2.59 and 
-4.98 x 10", respectively. A diagram of the perturbation velocities 
is presented in Figure 1. for a graphical representation of the 


solution determined in this Chapter. 
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V. THs VARIABLE COZFFICIZNT CASE 


8 The Solution of the Perturbation. 

with D determined and constant throughout the atmosphere for 
any particular mean latitude, the second order differential equation (28) 
is now solvable, By suitable manipulation, described in the Appendix, 
this eguation (28) is & case of the confluent hypergeometric differ- 


ential eguation 
pL“) + Ee-pl Lip - o. Lep) = © 


» 2 | 一 xz]. 


Equation (55) is satisfied by the functions 


(55) 





Y = F(a, bp) 21+ Sor 5 Gd 2 


٢ yo F(a-b+l,2-b, p) = 


t+ A-btı + L (a-bi I (a - e * a 6 
229 7 CA-b)(S-b) ? - 


The confluent hypergeometric equation then has the general 


solution 


HD = CY +CH, (57) 


(27) 
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To test the convergence of Y -FCa,b,p) , compare to the 
series E* , where X=-%% lol . It can be shown that, termwise, 
ዛ < <* and orders higher than the second can be dropped with 


accuracy to the first five significant figures. Then 


= a. - adat) 
٨ [+ b p + E Teo" vr 


where afo = ሥሪ፡-| ， 
Testing the convergence of 
سر‎ Y, = ጃሻር6-»ዛእወ-ኳ 9] 


after the first 12 terms, the series is more rapidly convergent than 
where Xa - af bl In the range of values used in this study, however, 


there is negligible error, if series is shortened to 


መሠሪያ FE, iş bel pL (a-b+)(a-b+ pa (59) 
* ^ (a-b)(3-v) 


The general solution in the troposphere from equation (57) is 


HG) = Ci C It apt LALAH) or) + 


Clo tt) 


Crp)? (14 O=b+l 4 „1 Q-b+iXQ-b+2), 2] , (60) 
3 m 2-b 7 2 (a - b) 3-b) D 


(28) 
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— Aoplication of the Tropopause Boundary Conditions. 

After the manner of evalusting coefficients in Chapter IV, 
equating 6 er and horizontal velocity across the tropopause, 
two resultant equations can be solved simultaneously for C, and 0 . 
D. The Complete Solution. 


Jith C4 and C3 determined, the perturbations in the troposphere 


are defined as 


u, CGP, HG LAD "red p'Q-9 9 |) sm mace?) 


— m (C, 0 j. Cs LA Cos m Cx-ct), 


K 


5 


a E ditt dei Sir? 4 m(x-< t), 


(61) 


g = م‎ [l-«z] የና + Cs f(y) 4 + 
Cp) ^p (l-b) Vi) پل‎ ረ Y +C, OI Sın mü-cE) 


= Dm (2m - cl 1 
? ex K | 


(29) 
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In (61) the primes denote differentiation with respect to z. 6 
values of parameters to be inserted in (61) are 


D Bo 75 n = اي1‎ x 10-8 enz, o = l9 x 1074 emi, 


K = 1.064 x 1979 5”? ) so that p can be evaluated for any value of z. 


In particular, at the tropopause 2 has the value Y = -.0176. Like- 
wise y = 415 x 1079, and is constant throughout the troposphere. 
The complete solution includes the solution above the tropopause 


which remains precisely as developed in Chapter IV-2. 


(30) 
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VI. DISCUSSION OF RESULTS 


It would appear that the variations of the value of no determined 
by the discriminant of the stratospheric auxiliary equation (32) plays 
the most important role in this entire investigation. By variation of 
this parameter alone, the wave velocity, the stationary wave length 
and the ratio of horizontal to the vertical maximum velocity, can be 
varied. within the allowable range of this parameter, wave velocities 
determined by Rossby are included. 

Gravity and sound waves are also included in these solutions for 
appropriate values of n. The values of n correspond to the inequality 
He. | < ይባ h, :4mh Z 25 , with values of the wave-length, I, 
suitably chosen to be representative values of wave-length for sound 
and gravity waves, respectively. It will be remembered that Charney 
15] ; in dealing with the problem of the unstable baroclinic wave, 
first eliminated these short stable waves, which in the present paper 

is equivalent to restricting the condition to D « مل‎ 

A further possible investigation is suggested by this study. By 
limiting the lateral extent of the disturbance, larger values of the 
velocity components are consistent with the small values of vertical 
divergence used in this study, the values of 4ይ 
would then partially compensate one another, giving the proper (small) 


magnitude to Y, Y = = >. š 


ou) 
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APPEWDIX I. 


The transformation of equation (28) to the form of a confluent 
hypergeometric differential equation follows. 


Zguation (28) is restated here. 
(i- x2) H@) -K HG + [^ 80 xz) | u) ل‎ (28) 
Divide by (A+B): 


l~- az SMS i Baz Çi O, (28.2) 
p» hs Ks nto (r- 2 H 


Divide byr/t. then multiply b A+ B 
ivide y(t Baz A+B) Z ره‎ en multiply by (A+ B) 


I- q2 K 
| — Baz Hiz) > 32) HC) + (A+B) HC) -= o (28.5) 
A+B A+B 


Momentarily, consider the first term only 


521 Baz. መ - A+B +] 
F Ozz fF AtB B 
AtB کی‎ Baz 
"FO 
AtB ER: 5 (2886) 
心 O(- Baz \ j 
AtB 


025) 





Expending, 


Baz gel É Baz ; (28.4) 
AtB ax | 


if it is assumed that higher order terms negligible. 


Souation (28.b) together with (28.c,d.) yield 


AtB LA Së 2 7 | j — 
4 (Is Bez) — era pe | ws A+B| He) =0. call) 


Collecting terms, (28.e) is 


[253 ሀ HE E. dd + A H&) + 14+ 可 H@) =6. er) 


Let 


5 = l- Aaz AZ E _ 
AtB? MB ` I- Ë KSE ` 


then 


| ተር6፳ .. -8 
A+B +> > ) 


and equation (20.f) is 


$ HG, - K(IF&-$) Hé) 4(A*B) Hso, (28.8) 


(36) 





where 





2 3 2 2 . ሪ 2 
si As) = J(9 44. 


Equation (28.2), in terms of $ , becomes 


mo Ich )شر‎ 45)2t6 +(A+B) JG) = o. — 
Dividing (28.h) by (ax کل‎ "٢ 


6 [A+B] 


By letting re - K B(AB)/A*4 
A 


equation (28.i) becomes 


777 و دري 7 *(4+8) م روغ ~ r (AtB‏ + (ع) ل ؟ 
f r (At 6) SC) + 72 jop-*.‏ 


(27) 
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Then if WË = P, 
Teen. — J (p. LEI = L), 
64ل‎ = 


1 
sr 
E 
4 
ን 
rm 
S 


3 03) 


eguation (28.3) is 


M 


“ብዬ p)" vL», 


Al d 2 4 P = ° 
BrLy) rr(Azl - 5 كا(‎ 4 Brey L(p)=0. G9 


Divide equation (2&.k) by ም ana substitute the equality p = Fa , 6 


the equation becomes 


۹ ار ے tr D -ol L oe‏ غو 
B | h 2‏ | 
If bz 一 (4+B)AB , Ue AB (any :‏ 


equation (28.1) is identically eoustion (55), which is the form of 


the confluent hypergeometric differential equation. 


(38) 
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